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^ '. Abstract 
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' The structure of the Euler-Lagrange equations for a general Lagrangian theory 

■ (e.g. singular, with higher derivatives) is studied. For these equations we present 

a reduction procedure to the so-called canonical form. In the canonical form the 
equations are solved with respect to highest-order derivatives of nongauge coordi- 
Q"^ ' nates, whereas gauge coordinates and their derivatives enter in the right hand sides 

OO . of the equations as arbitrary functions of time. The reduction procedure reveals 

I constraints in the Lagrangian formulation of singular systems and, in that respect, 

is similar to the Dirac procedure in the Hamiltonian formulation. Moreover, the 
^ . reduction procedure allows one to reveal the gauge identities between the Euler- 

Lagrange equations. Thus, a constructive way of finding all the gauge generators 
within the Lagrangian formulation is presented. At the same time, it is proven that 
for local theories all the gauge generators are local in time operators. 



> 
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1 Introduction 



X 

^ • At present increasingly complicated gauge models are used in field and string theory. Gen- 

erally a comprehensive analysis of their structure is not a simple task. In the Lagrangian 
formulation the problem includes, in particular, finding generators of gauge symmetries 
and their algebra, revealing the hidden structure of the equations of motion and so on. One 
ought to say that in the Hamiltonian formulation there exists a relatively well-developed 
scheme of constraint finding (Dirac procedure 0]) and reorganization ^ ^, The 
constraint structure can be, in principle, related to the symmetry properties of the initial 
gauge theory in the Lagrangian formulation [^]. However, in the general case, this relation 
cannot be considered as a constructive method to study the Lagrangian symmetries (it is 
indirect and complicated). Moreover, the modern tendency is to avoid the non-covariant 
hamiltonization step and to use the Lagrangian quantization [^] for constructing quantum 
theory. Such an approach incorporates all the Lagrangian structures (in particular, the 
total gauge algebra). That is why it seems important to develop a reduction procedure 
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within the Lagrangian formulation - in a sense similar to the Dirac procedure in the 
Hamiltonian formulation - that may allow one in a constructive manner to reveal the 
hidden structure of the Euler-Lagrange equations (ELE) of motion and to find all the 
gauge identities and therefore the generators of all the gauge transformations. An idea of 
such a procedure was first mentioned in publications of the authors (D.G and I.T) 0, || 
(see also Appendix C in [^), but was not appropriately elaborated and some important 
points where not revealed. 

In the present paper we return to this idea studying the structure of the ELE for a 
general Lagrangian theory (singular, with higher derivatives, and with external fields). In 
Sect. II we introduce some notation and definitions. In Sect. Ill, we reduce the ELE of 
nonsingular theories to the so called canonical form (in the canonical form the equations 
are solved with respect to highest-order derivatives of nongauge coordinates, whereas, 
gauge coordinates and their derivatives enter in the right hand sides of the equations as 
arbitrary functions of time, see below). In Sect. IV we formulate the reduction procedure 
for the singular case. In a sense, the reduction procedure reveals constraints in the 
Lagrangian formulation of singular systems and, in that respect, is similar to the Dirac 
procedure in the Hamiltonian formulation. In Sect. V we demonstrate how the reduction 
procedure reveals the gauge identities between the ELE. Thus, a constructive way of 
finding all the gauge generators within the Lagrangian formulation is presented. At the 
same time it is proven that for local theories all the gauge generators are local in time 
operators. In the Appendix we collect some Lemmas being useful for our consideration. 



2 General ELE 



2.1 Notation, definitions, and conventions 

We consider a system with finite degrees of freedom (classical mechanics). These degrees 
of freedom are described by the generalized coordinates g", a = 1, ...,n, which depend on 
the time t. The following notation is used: 

g'^ = ^,g" = ^,---, °rg'^W = ^, / = 0,1,..., (g»ra = g"). (1) 

The coordinates g° = are called sometimes velocities of zeroth order; the velocities 
ga _ g,a[i] g^j^g called velocities of the first order; the accelerations = are called 
velocities of second order, and so on. The space of all the velocities is often called the jet 
space, see 0. 

As local functions (LP) we call those functions that are defined on the jet space and 
depend on the velocities g"''' up to some finite orders A^a > (/ < No). Further, we call 
Na the order of the coordinate in the LP. Por the LP we use the following notation^: 

F {q\ (f, if, ...) = F (g*^™, g^'^U"'^', •••) = F (gW) , g^ = (g'^W, < / < A^J, 

or sometimes : F (g^) = F {■ ■ ■ q"^^^^) . (2) 

In the latter form, we indicate only the highest-order derivatives in the arguments of the 
LP. 



""^The functions F may depend on time explicitly, however, we do not include t in the arguments of 
the functions. 
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The following notation is often used: [a] is the number of the indices a, namely, if 
a = 1, ...n, then [a] = n. Similarly, suppose Fa (rj) , a = 1, ...,n are some functions, then 
[F] is the number of these functions, [F] = n, etc. . However differently, writing we 
denote by [/] the order of the time derivatives, see (|1|). 

On the jet space, we define local operators (LO) to be matrix operators U of the form 

(3) 



fc=0 

,k 



where K is a finite number and u^a are some LF. The LO act on columns of LF /„ 
producing columns of LF Fa = UAafa as well. We define the transposed operator to U as 

r).,-i:B)'''^.-i:H)'E(:)''S'(|)"- w 

A;=0 ^ ^ fe=0 1=0 ^ ^ ^ 

Then the following relation holds true 



FAUAafa 



Un Fa 

aA 



fa + ^Q. (5) 



dt 

where Fa , fa , and Q are LF. The LO Uab is symmetric (+) or skewsymmetric (— ) 
whenever the relation (^U^^ = ±Uab holds true. 

Suppose a set of LF Fa ■ ■ g^t^" ] j , or a set of equations Fa ■ ■ q^^^" ] j = , be 

given. In the general case the orders A*"^ of the coordinates q'^ in the functions Fa (in the 
equations Fa = 0) are different, i.e. these orders depend both on a and A. The number 
Afa = max^ is called the order of the coordinate in the set of the functions Fa (in 
the set of the equations Fa = 0). 

Whenever, a subset Fa' = 0, A' C A has orders Af^ of the coordinates less than the 
corresponding orders of the complete set, namely, Va : < Afa , we call this subset the 
constraint equations. 

Generally two sets of equations. Fa (g^'') = and fa (g^'^) = are equivalent whenever 
they have the same set of solutions. In what follows we denote this fact as: F = <^=^ 

/ = o. 

Suppose that two sets of LF Fa (g''^) and Xa (q'^'O > [F] = [x] > are related by some LO, 

F = Ux, X = VF, UV=1. (6) 

Then we call such functions equivalent and denote this fact as: F x ■ Obviously, 

F~x^i^ = O^X = 0. (7) 

If d^) holds true, we will call the equations F4 = and = strong equivalent. 
In what follows we often meet the case where 

Xa=(o^),^=(«,G');VG:Og = 0. (8) 

Here the equivalence F ~ x implies the equivalence of the equations F = and / = 
and the existence of the identities VqaFa = 0. Namely, 
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2.2 ELE 



Below we restrict our consideration to the Lagrange functions L that are LF on the jet 
space, and depend on some external coordinates (fields) u'^ (we call the coordinates 
external ones in contrast to the coordinates , which we call inner coordinates) which 
are some given functions of time. Thus, 

L = L (■ , a = l,...,ri, A^^ > 0. (10) 

The orders A^^ of the inner coordinates q'^ in the Lagrange function will be called further 
the proper orders of the coordinates. Coordinates with the proper orders A^^ = 0, we 
call the degenerate coordinates [|10]- 

Equations of motion of a Lagrangian theory (the ELE) follow from the action principle 
5S = 0, S = J Ldt , and have the form (merely the inner coordinates have to be varied): 

l^ = Ef-TV#m=0' a=l,...,n. (11) 

Following |jlO| , we classify the Lagrangian theories as nonsingular (M 7^ 0) and singular 
(M = 0) ones by the help of the generalized Hessian M = det ||Mafc||, where 

^"'^ " d^fmo^^] ^^^^ 

is the generalized Hessian matrix. 

In what follows the ELE of a nonsingular (singular) theory will be called the nonsingular 
(singular) ELE. 

Sometimes, it is convenient to enumerate the inner coordinates and organize them 
into groups such that = (g"o, ...,q°-') , where are groups of indices that enumerate 
coordinates having the same proper orders, Na, = . Besides, we organize these groups 
such that rii > rij-i ■ ■ ■ > uq = (max Na = Naj = nj , and are the degenerate 
coordinates, Na^ = uq = 0). Thus, 

a= {ttk, k = 0, !,...,!) , [a] = ^[ai], [ai] > , rij > rij^i ■ ■ ■ > Uq = . (13) 

i 

Taking into account the notation (0), we may write the Lagrange function and the ELE 
as: 

L = L (■■■g'^'=["*];M^) , A; = 0, 1, / ; (14) 
Fa, (■■•g^[^''+"*]; =0, (15) 

Here Ma, b is the generalized Hessian matrix and Ka, and M^y are some LF of the indicated 
arguments. 

Consider the orders of the inner coordinates in the complete set of the ELE. These 
orders are A^i = Na + rij . One can see that these orders are, in fact, defined by a subset 
of ( P^^D with = / . In any subset of the equations (|T^) with k < I the orders of the 
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coordinates are less than in the complete set. Then according to the above definition, all 
the ELE with k < I are constraints. The set (|T5|) has the following specific structure: In 
each equation of the complete set the order of a coordinate is the sum of the proper 
order Na and of the order Uk ■ The latter is the same for all the coordinates and depends 
only on the number of the equation. 

2.3 Canonical form 

Let a set of equations 

=0, (17) 

be given. Suppose that these equations can be transformed to the following equivalent 
form: 

g-M = (. . . . ^ qa ^ ^^a^^g^ ,a={a,g), la<Ua. (18) 

The equation (|T^) present the canonical form of the initial set ([T7|). In the canonical 
form the equations are solved with respect to the highest-order time derivatives g^I'^l of 
the coordinates q". The coordinates q^ (if they exist) and their derivatives g^t'^l enter 
into the set (0) as arbitrary functions of time. In fact, there are no equations for these 
coordinates. In what follows we call these coordinates the gauge coordinates whereas 
we call the nongauge coordinates. The orders of the coordinates in the canonical forms 
may be less than those in the initial set. In the general case, one and the same set of 
equations can have different canonical forms. Generally there are many canonical form of 
the same set of equations. 

Below, we are going to formulate a general procedure of reducing the ELE to the 
canonical form (in what follows it is called the reduction procedure). Our consideration is 
always local in a vicinity of a given consideration point ^q^'^ (in the jet space), which is on 
shell w.r.t. the ELE. We consider theories and coordinates where the consideration point 
could be selected as zero point. Thus, we suppose that the zero point is on shell. Further 
we always suppose that the ranks of the encountered Jacobi matrice^ are constant in a 
vicinity of the consideration point. Such suppositions we call "suppositions of the ranks". 
Saying that some suppositions hold true in the consideration point, we always suppose that 
they hold true in a vicinity of the consideration point. In course of the reduction procedure 
we perform several typical transformations with LP or with the corresponding equations. 
Each of such transformations lead to equivalent sets of equations or to equivalent sets of 
LP (definitions of such equivalences are given above). The proof of these equivalences is 
based on two Lemmas which are presented in the Appendix. Any statement of the form 
"the following equivalence holds true" can be easily justified by these Lemmas. 

3 Canonical form of nonsingular ELE 
3.1 A particular case 

Consider theories without external coordinates and with only two different proper orders 
of the inner coordinates. In such a case all the indices a can be divided into two groups: 

retangular matrix with elements dAajdx^ is often denoted as dA/dx and called the Jacobi matrix. 
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a = (ai, aa) , such that A^'^^ = na > A^^i = , L = L [■ ■ ■ g^^al^^l, ■ ■ ■ g^il'^il) . Consider first 
the case > 0. Then Eqs. (|I5[) can be written as: 

Fa, = M,,,g''[^"+"^1 + Ka, (■ ■ ■ gMA^.+n2-l]) ^ Q , (19) 
i^ai = M,, + (. . . gMA^.+"i-l]) = . (20) 



The equations ( pOf ) are constraints. Consider the set 



(21) 



obtained from the constraints after being ria — Ui times differentiated with respect to the 
time t. Since M 7^ 0, the rectangular matrix Ma^ a has a maximal rank, that is why there 
exists another division of the indices: 



(a|i) , [a|i] = [tti] , z = 1, 2 , det Ma^b^^ ^ . 



Remark that 



fli = (aj|i,ai|2) , a\i = (ai|j,a2|i) , [0112] = [a2|i] ■ 
The set (BTl) can be solved with respect to the derivatives follows: 



(22) 
(23) 

(24) 



Here the matrices Mi and M3 are defined by the following block representation of the 
matrix M: 



M, 



det M„jb ^ ^ det Mi ^ 



Excluding the derivatives g^iH^"!!"*""^] from Eqs. (|T9|) by the help of (p4|) , we get the 
equations 

(M5),,,,^g'''^KH (...gMA^^+--y) =0. 

Taking into account an useful relation 



(25) 



^^^^ = ^"*^Mi M3 J=^"'^Mi Ma 
= det Ml det (M4 - MgMf ^Mg) , 

which is related to the Gaussian reduction of matrices 0], we get: 



(26) 



det M ^ I ^ _^ Ms = M4 - MaMf ^M; 
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det Ml ^ 

Therefore, (E^) can be solved with respect to the highest-order derivatives ^"^12 [^"12+" 



(27) 



as: 



q-l2[Na,,+n2\ ^ _(^M^iyi2a2 ;^^^(...gM^.+n2-i])_(M2Mr^);;^i^(^)(---g''[^''+'^^-i]) 



b|2[A'6|2+"2-l] _ . . „^|l[^b|i+"2-l] 



= I . . . g 'n "|2 — J , ■ ■ - g 



(28) 



Thus, we get a set 

F,^ = (MO^^^i^g'^i^KiH =0, (30) 

which is strong equivalent to the initial ELE by virtue of the Lemma 1. 

Due to the condition det Mi ^ 0, the equations ( pOD can be solved with respect to 

and we obtain: 



(31) 



Eqs. (^) and (^) are not of canonical form since the functions (/^^i^ contain derivatives 



Nb.^ +ni-l 



. Now we exclude all the 



^''|i[-'^f'ii+"2 exceeding the "allowed" order 
surplus derivatives g^i^f^"!!"*""^] ^ ...^ g"ii[^"|i"''"'^ from the right hand side of (|29| ) by the 



help of (plD and corresponding derivatives of it. To this end we need to differentiate (|3l| ) 
not more than n2 — ni — 1 times. Finally, we obtain the following strong equivalent form 
(the equivalence is justified by the Lemma 1) of the ELE: 



a|i[Af,|^+ni] ^ ja|i ( . . . g''|2 [^''|2 . . . q^\^[^t'^i+'^l~l] 



(32) 



It is just the canonical form. Taking into account the division of the indices w.r.t. proper 
orders of the coordinates, one gets: 

^a2|2[2n2] _ ja2|2 ^. . . gb2|2[2n2-l] . . . g6l|2["l+"2-l] . . . gfe2|2["2+ni-l] . . . g?*!!! P"! -1] ^ 
g«i|2["i+"2] _ jai|2 ^. . . g^'2|2[2n2-l] . . . qbi\2lni+n2-l] . . . gf)2|i ["2+ni -1] . . . gfei|i [2ni -1] ~J 
^a2|i[ri2+ni] _ ya2|i ^. . . g^'2|2 ["i+n2] . . . g''i|2[2ni] _ _ _ g62|i["2+ni-l] . . . qh\i['^rLi~l]^ 
gai|i[2ni] _ jai|i ^. . . ^^212 [»^i+»^2] . . . gbi|2[2"i] . . . ^^211 [n2+rii-l] _ _ _ [2rii-l] ^ j^gg-j 

Remark that the number of the initial data is equal to 2 Na . Indeed, 

[0212] (^2 + + [ai\2] (ni + 722) + [a2|i] (n2 + ni) + [ai|i] (ni + ni) 
= 2[a2] n2 + 2[ai]ni = 2j2Na. 



One ought to mention that the canonical form (^) was obtained in [jT2|. However, the 
procedure that was used for that purpose did not provide the proof of the equivalence 
between the initial ELE and the form (B^). 
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Suppose now that the Lagrange function contains degenerate coordinates a = 
(ao, ai) . Thus, L = L (g"", ■ ■ ■ and the ELE read: 

Fa, = Ma, + Ka, (■ ■ ■ g^I^^+^i-H) = , (34) 

^ :M.„(...gM^^l)=0. (35) 



ao 



dq" 



Despite these equations are formally different from the above case, the whole procedure 
of reductions goes through without any essential change. In fact, differentiating Eq. 
111 times, one obtains 

M„„„g'^[^»+"^1 + Ki'J (■ ■ ■ = 0, (36) 

and all the previous steps may be done as before. Namely, one obtains 

g^l^KlH =^'^11 /^...g''|lKl+"l-l],...gVK+"«-^]^ , (37) 



and, since det ||(Mi)aQ(jnj|| 7^ 0, Eqs. ( P^D can be solved with respect to the variable 
g^'^t^^ioJ as follows: 



Finally, after eliminating the "bad" derivatives in the right hand side of (|37|) for g'*'^ 

one ends up again with Eqs. ( ^3]) but now with ^2 — > ni, ni — >■ (by convention: 

qbnii-n = 0). 



3.2 General nonsingular ELE 

Consider the general nonsingular ELE. Here the Lagrange function may contain some 
degenerate inner coordinates, higher derivatives of some inner coordinates, and, moreover, 
may depend on some external coordinates, L = L (■ ■ ■ 

^a[JV„].^/.j , Na>0. Thus, we are 

going to deal with the nonsingular ELE of the form (p^). Our aim is to present these 
equations in an equivalent canonical form. 

Theorem 1: The nonsingular ELE ([TBI) can be transformed to the 
foUowing equivalent canonical form: 

/>A;+>A; + 1, k > > 0, i, j, k = 0,1, I , (38) 

where the indices of the coordinates are divided into groups as follows: a = 
(tti) is the division of the indices w.r.t. the proper orders of the coordinates, 
and besides 

fli = {ai\k, hk = 0,1,...,/) , [ai\k\ > 0, ^ [0^1^] = ^ [ak\i] = [ai] = [a\i] . 

k k 
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Moreover, the equivalence F ~ / between the corresponding LF holds true. 
That implies 

Fa = Uabf , f = V'^Fa , UabV"' = 5^ , 

where U and V are LO. Besides, that implies the strong equivalence between the 
ELE and their canonical form (p^). 

The proof of the Theorem 1 may be considered, in fact, as the general reduction 
procedure to the canonical form for the nonsingular ELE. 

It is reasonable to divide the reduction procedure into two parts. These parts may be 
called conditionally "the preliminary resolution", and "the subordination procedure". 

Preliminary resolution 

Let us introduce the notation a = (a, a/) , a = (a^ , A; = 0, 1, ...,/— 1) , Na < nj , 
such that the ELE read: 

Faj (■ ■ ■ g''[^''+"^]; ■ ■ ■ = Maj + Ka, (■ ■ ■ g^I^^+^^-i]; . . . m'^I"^]) = , (39) 

F^('...gf'K.+^J;...^4^J) =0. (40) 



Recall that the equations ( ^OD can be considered as constraints. 

The first step of the procedure is the following: We consider the consistency condi- 
tions of the constraints. Namely, we consider the equations that are obtained from the 
constraints by differentiating them nj — Ua times, 

^h-^J = M,,g^[^''+"^] + ir« (■ • ■ gMiV.+n.-l]. . . . ^M".]) = . (41) 

Here K^P are some LF of the indicated arguments. Remark that the orders of all the 

coordinates in the set ( |41| ) coincide with the ones in the complete set. For M 7^ 0, the 
matrix 

Qgb[N,+ni] = Qqa[Na]Qqb[N,] = (^2) 

is invertible. At the same time, the rectangular matrix Maa has the maximal rank [a]. 
Therefore, there exists a division of the indices a such that: 

a = (a, a\j) , [a] = [a] , [a\i] = [a/] , det Maa 7^ . (43) 



Thus, the division (|TB|) of the indices a w.r.t. the coordinate proper orders becomes more 
detailed, 

tti = (oi, ai|/) , a = (ai) , a\i = (0^1/) , [ai\i] > , 

= hi] = [«^] ' Yl [^*] = [^] = [-] ■ 

i i 

Due to (^3]), the set (^Tj) can be solved with respect to the derivatives 



, (44) 
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where 



M, 



Excluding the derivatives g'^[^°+"-r] f^om Eqs. (^) by the help of (0), we get the set: 

Ms = M4 - MsMf ^Ms , det Ms ^ , (45) 

where k!^j^ are some LF of the indicated arguments. The set (|45D can be solved with 
respect to its highest-order derivatives 



(46) 

where (p°'^' are some LF. Thus, after the first step we get a set of equations 

^ _^a^, ..gf'iiK,+«i-i]^...g5[W5+«/-i]....^/'M^ =0, (47) 

=0, a = (afc, A; = 0,1,...,/ -1) , iV, < iV, , (48) 

which are strong equivalent to the initial ELE by virtue of the Lemma 1 from the 
Appendix. 

At the second step we turn to the subset (^Sj) . We remark that this subset has the 
same structure as the complete initial set of the ELE if one considers the coordinates 
as inner ones and the variables as external ones. Namely, let us denote 

F^(^... gSh+^J ; . . . „Mi [N^ ) = (■ ■ ■ g^K+^J ; . . . 



Then the set (ff3) can be written as: 



>^(^...g"^K5+A'J;...„MiK]j =0, a = (a,. A; = 0,1,...,/ -1) , Na< Nj , (49) 



where 



Mao (■ ■ ■ g^f^^^; M^i) = dL/dq"" 



Here are the inner coordinates, and m^i are the external coordinates. The order of the 
set ( ^91) is 2?2/_i. Furthermore, by virtue of (^3]), the matrix 



nr_i -Na 



1 r-r-i 
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is invertible. Thus, the structure (|T5|,16) is repeated completely. 

At the same time, the number of the inner variables, the number of the equations, and 
the order of the set (|49|) are less than those of the initial set of the ELE ([T5| , p^ . 



Now, we apply the same procedure as in the first step to the reduced set p9|). That 
will be the second step of the reduction procedure. It will produce equations of similar 
structure with less inner variables and of lower order. After the last (/ + l)-th step the 
ELE (|15]) may be written in the following strong equivalent form: 



J > A;+ > A; + 1, > > , J > z,j > 0, (51) 

where 0"'''' are some LF of the indicated arguments (the arguments ■ ■ ■ result 
from those coordinates that intermediately have been considered as external ones), and 
the indices are divided into the following groups: 

ai = {ai\k), [ai\k]>0, ^ [ai\k] = ^ [ak\i] = [ai] , i, k = 0,1, I . 

k k 

The set (^) is still not the canonical form of the ELE. The reason is that the right 
hand sides of the set contain derivatives of orders that may exceed the orders rii + Uk of 
the (highest) derivatives appearing on the left hand side of the set. We recall 

that by the definition in the canonical form there is a subordination of derivative orders, 
namely, the orders of all the derivatives in the right hand sides have to be less than the 
ones on the left hand side. Explicitly, this subordination would require that the following 
inequalities should hold: 

rij + rik^ > Uj + rik , 
rij + nk_ > rij + Uk - I , 

which, because of the inequalities ni > > ■ ■ ■ ni > no , is true for the first line and 
the case k- = k of the second line, and it is definitely not true for the cases k_ < k. 
Arranging the equations ( p^ (for fixed value of i) in descending order w.r.t. k, and the 
arguments in the functions (f (for fixed value of j) also in descending order w.r.t. the 
value of k^ and k_, we get, when disregarding the common value rij, a quadratic matrix 
whose main diagonal (i.e. elements with k = k_) contains the entries — 1, whereas the 
entries to the left of that diagonal are equal to n^, and to the right of that diagonal are 
equal to nfc — 1. Therefore, below the main diagonal occur "good" derivatives, and above 
it occur "bad" derivatives not obeying the subordination requirement. 

Subordination procedure 

One can see that these "bad" derivatives can be excluded from the right hand sides by 
the help of corresponding "lower" equations of the set and their differential consequences 
(compare Eqs. ( p9D and (0) for the simple case 1 = 2). In what follows we call such an 
exclusion the subordination procedure. 

In order to be more definite let us write down two arbitrary lines, i > k, of the right 
hand sides of the set of equations (0) (for the highest derivatives only): 
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Obviously, because > all the derivatives of the equation for with 
k > £_ > are "bad" with respect to the derivatives g'^iifc["-»+"'=l (remind £ > £_ > 0). 
However, these "bad" derivatives can be eliminated by the equations for the latter ones, 
q"-i\klni+nk] ^ and their differential consequences up to the order n£ — Uk — 1. Thereby, the 
function 0"''^ changes into some new function One can see that doing this we do 
not change the highest orders of derivatives of the other coordinates, both proper and 
external ones, in the right hand side of the equation for (Remind, that the 

derivatives of the external coordinates are and m'^'"*!, respectively.) 

This subordination procedure, starting with i = I may be done for any k < I, thereby 
"cleaning" every entry on the right hand side of equations for Namely, the 

highest orders of derivatives on the r.h.s. become g^'ji*- with / > A;_ > (for 
the case i = I no k+ appears). Then the procedure will be applied to the equations for 
^ai|j_i[jii+n7_i]^ and so forth, up to g"i|o["i+"o]^ where nothing is to be changed. 

After having eliminated all the "bad" derivatives, we transformed the set (|5TD, and 
therefore the initial ELE, to the following strong equivalent (the equivalence is justified 
by the Lemma 1) canonical form 

I >k+>k + l, k>k_>0, i,j,k = 0,1,...,/, 

where y^"'!*^ are some LP of the indicated arguments. This proves the Theorem 1. 
We see that there are no gauge coordinates in the nonsingular ELE. 
The number of the initial data is equal to 2 J2a ■ Indeed, 

i,k i \ k / k \ i / a 

One ought to remark that in the general case there exist many different canonical forms 
of the nonsingular ELE. This uncertainty is related to the possibility of different choices 
of nonzero minors of a matrix with a given rank (different divisions of the indices a, in 
course of the reduction procedure). However, as it was demonstrated above, the number 
of the equations in the canonical form (which is equal to the number of the ELE in the 
nonsingular case) and the number of the initial data is the same for all possible canonical 
forms. 



4 Canonical form of singular ELE 

Studying the canonical form of nonsingular ELE, we have demonstrated that the equations 
in the canonical form are solved with respect to the highest-order derivatives g'^i|fe["-»+"-'=l^ 
where Ui are the proper orders of the coordinates q"'' . However, considering specific 
examples, one can see that this is not always true for singular ELE. Namely, in the 
canonical form of the latter case, the highest orders of the derivatives g"*''' may take on 
all the values from zero to + / . The reduction procedure to the canonical form for the 
general singular ELE is considered below. In the singular case, already after the first step 



of the reduction procedure, the ELE cease to have their initial specific structure (15,T6| 



Namely, the simple structure of terms with highest-order derivatives in the equations may 



12 



be lost. That is why in the singular case it is more convenient to formulate the reduction 
procedure for a more general set of ordinary differential equations, which contains the 
ELE as a particular case. Namely, further we are going to consider a set of the formQ: 

^^^(...^^4.+/.]) ^0; 2 = 0,1,...,/, /i = 0,..., J. (52) 

Here Fa^ (■ ■ ■ are some LP. Via and are denoted sets of indices, [oj] > 



0, > , and the complete set of the inner coordinates in Eqs. (|5^) is g° 

g"^) , a = (oj , i = 0,1,...,/). The indices A = (A^) enumerate the equations. 
In the general case the number of the indices A (the number of all the equations) is not 
equal to the number of the indices a (the number of the coordinates). The division of the 
indices A into the groups is not related to the division of the indices a into the groups. 
The orders of the coordinates q"'^ in the complete set (|5^) are: A/^. = i + J . In fact. 



these orders are defined by a subset of (^) with /i = J. In all the other equations with 
fi < J the coordinates q""' have the orders less than i + J. Thus, the latter equations are 
constraints. 

Similar to the ELE ([T^p, the set ( [5^ ) has the following specific structure: In each 
equation of the set the order of a coordinate is the sum of the proper order i and of 
the order /i. The latter is the same for all the coordinates and is related to the number 
of the equation in the set. 

Below we consider the reduction procedure to the canonical form for the equations 
(0). In fact, this reduction procedure is formulated in the proof of the Theorem 2 given 
below. The Theorem 2 holds true under certain suppositions of the structure of the 
functions Fa,, ■ These suppositions are formulated as suppositions of the ranks of some 
Jacobi matrices involving the functions Fa^, ■ First of all, the complete matrix 

has to have a constant rank in a vicinity of the consideration point (one can see that the 
matrix MA^at coincides with generalized Hessian matrix if the set (0) is the Lagrangian 
one). 

Theorem 2: Under certain suppositions of the ranks, the equations 
( |521) can be transformed to the following equivalent canonical form: 

j = 0, 1, ...,/, cr = -/, ...,J, -/<(T_<a,cr + l<(T+<J + l, (54) 
where all the indices a are divided into groups as follows: 

ai={ai\a), [ai\a]>0, a = -/,..., J + 1, ([ai|^] = if i + a < O) , (55) 

and it is thought that negative powers of the time derivatives do not exist, 
that is: = for p < 0. 

■^We do not indicate here possible external coordinates. 
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Moreover, the following equivalence between the corresponding LF holds 
true: 



Fa~Fa= ( ^J^'' ) ' (^ik.G) , ^ = 0,1,..,/, a = -/,..., J, 

Og = VG , [G] = [A] - [a] + [a,\j+,] . (56) 

i 

That implies 

Fa = U^Fb , Fb = V^Fa , U^V^ = 5^ , (57) 

where U and are LO. 
Let us make some comments to the Theorem 2. The canonical form ( |54l ) of the singular 
ELE differs from that ( ^81) of the nonsingular ELE. As was demonstrated in the previous 
Sect., in the latter case the spectrum of the orders of the variables in the canonical 
form extends from i + yUmin to i + J . In the singular case, we have to admit (and one 
can see this on specific examples) the spectrum extends from to i + J . Under such a 



supposition we can justify by the induction the structure (|5^) of the canonical form. One 
can see from (^3]) that each group of the indices Oj is divided in subgroups — > ai\„ , 
a = —I, J + 1. In the canonical form the singular ELE are solved with respect to the 
highest-order derivatives g"'i<'[*+'^l , a = —I, J, ( [aj|o-] = for i + a < 0). There are 
no equations for the coordinates g"»iJ+i . These coordinates enter the set (|54l ) as arbitrary 
functions of time. They are gauge coordinates according to the general definition. As in 
the nonsingular case, it is supposed that no coordinate in the function y^'^'i'^ has an 
order greater than k + a (the proper order plus a) . Besides, the order of the coordinates 
q''k\a_ function 99"*!'^ has to be less than k + a- . 

We are going to prove the Theorem 2 by induction w.r.t. M = I + J . To this end, 
we consider first equations of lower orders, then we use an induction to prove the general 
case. 

4.1 Equations of lower orders 

Remark that the case A/" = implies I = J = and the set (0) is reduced to form 

Fa (q) = 0, q= (g'^) . (58) 

Here the Theorem 2 holds true by virtue of the Lemma 3 from the Appendix. 

Let J\f = 1. That implies either 1 = 1, J = OorJ = 0, J = 1 . Consider, for example, 
the first case. Here (i = 0, 1, /x = 0) and the set (|52[) reads 



FA(g'^°,g'^\g'^0 = 0, N>0, [ao]>0. (59) 
In the case under consideration the supposition (|53| ) reads: 



rank — — 4t = f ■ (60) 

Then there exists a division of the indices: A = (74/1,74/2), = (04/1,01/2), [^/i] = 
[c^o/i] + ['^o/i] = such that 
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det 



/I 



7^0. 



Thus, we may solve the equations Fa,^ = with respect to g'^'/it 







(61) 



Then we exclude the arguments g'^viW from the functions F^^j by the help of ( |6lD 



By virtue of the Lemma 2 from the Appendix, the functions Fa^^ depend on only 
Thus, we have the equivalence^ 



(62) 



Now we suppose that the matrix dFA^^/dq"'^ has a constant rank. Therefore (see Lemma 
3) 



^ A/2 



, ai = (ai,ai) 



Let us exclude the arguments g-i, g-i from the functions -Fa/i, by the help of the equations 
Then the equivalence holds true: 



q^l =ip^ 



( 



F ^F 



ai ^aiN 



i^A/i {q''\q''\q 



a = (oq, Oi) , ai = {a^, ai) . 



(63) 



The set of functions F has the same structure as the initial set F. However, the number 

1 

of the nonzero functions F is less than the number of the functions F. Moreover, some of 

1 

the functions F depend linearly on a part of the variables. That is why the supposition 

of the type (^) for the functions F^^^ is reduced to the supposition about the rank of the 

1 

matrix dF a,^I d {q"'^ , (f^) . Accepting the latter supposition we apply the above reduction 

1 

procedure to the functions F and so on. After the i-ih stage we have the following 
equivalence: 



F ^ F 



FA/^iq^'^q^'^cf^) 

q^i _ i^a^ ^qai^ , a — (oq, cti) , cti — (flj, a^) . 



Here, and in what follows, we use Lemma 1 to justify the equivalence. 
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The procedure ends at a k-th stage when 



rank 



k 



/k 



Then there exists a division of the indices = (ai|o, o-gj , clq = (ao|0) So) > ['^i|o] + [c^oio] 
[A/fc] , such that 



det 



k 

dFA 



Ik 



7^0 



Denoting = ai|_i , G = Gfc, such that a = (ai|_i, ao|o, diio, , and a^, = (ago,agi) 
[G] = [A] — [a] + [flg], we get finally the equivalence: 

Og 



V 



(64) 



Here q"'^ = {q°-<^o , g'^si ) are gauge coordinates. Thus, the Theorem 2 holds true in this case. 
The case / = 0, J = 1 = 0, yU = 0, 1) corresponds to the equations of the form 



Fa, (g'^Sr) = 0, Fa, (g"i) = 0, 



(65) 



Such equations present a particular case ([ao] = 0) of the equations Fa = with the LF 
Fa defined in (|62D. The reduction procedure for the latter case was considered above. It 
leads to the following equivalence 



qO-\Q _ y9''|0(g"|i^ ^^9) 



(a[o,a[i,ag) , [G] = [A] - [a] + 



"gi ■ 



Here q""^ are the gauge coordinates. Thus, the Theorem holds true in this C3iSG clS well. 



4.2 Equations of arbitrary orders 

We have verified that the Theorem 2 holds true for A/" = 0, 1. Now we are going to prove 
the theorem for M = I + J = K (where K is some fixed number) supposing that the 
theorem holds true for any M < K . 
At the first step we consider the set 

^y^-M](...^a,[.+J])^0, ^ = 0,1,...,/, /i = 0,...,J, (66) 

which is obtained from the initial set (|5^) by substituting the constraints by the corre- 
sponding consistency conditions (conditions obtained from the constraints by J — /i 
time differentiations). According to the supposition (^3|), there exists a division of the 
indices A^ and as: A^ = [A^^/i, , = (0^/1,0^/2) , J^/, [^/^/i] = Ei = 
such that: 

^ . (67) 



det 



OF 



dq" 



i/i 



[i+j] 
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Thus, we may solve the equations ^'^^ = with respect to the derivatives . 
Namely, 



(68) 

bi/l[j+J] 



Now we pass from the functions Faj^^ ones Faj/^ excluding the arguments q^/^ 

from the former, 

Fa,,. . (69) 



- Faj/2 



lj/2 



The fact that the functions F^j^^ do not depend on both and g''>/2[*+-^l is based 

on the Lemma 2 from the Appendix. Thus, we have the equivalence (see Lemma 1 from 
the Appendix) 



V = 0,..., J- 1 



F\, 



A' ' 



(70) 



where 



Fa. , u = 0,...,J-2 



(71) 



Let us turn to the functions F^, . They have the same structure as in ( p2[ ) and 
correspond to the case M = I + J<K.\a accordance with the induction hypothesis, 
supposing, in particular, that the matrix 



M' 



A' a. 



9F'a, 



has a constant rank in the consideration point the following equivalence holds true: 



F'a, 



z,j = 0,l,...,/, [G'] = [A']-[a] + Y,[a,\j] , 

i 

cr = -/,..., J- 1, -/<cr_<(T, a + l<(T+<J-l. 
Taking into account (^), we obtain 



(72) 



v 



Og' 



i,j = 0,1,...,/ , a = -I,...,J -1, -/ < o-_ < a, a + l<a+<J-l. (73) 
Now we pass from the functions F^^,, to the ones F^,,, excluding the arguments g^'WlP'l 

J / 1 J / 1 

Pi > i + (J, (J = —I, J — 1 from the former. As a result, the following equivalence takes 
place: 



F ^ F 



Faj/^ (■ ■ ■ ■ ■ ■ qh\Ai+<^-A'^ 



f>jW_ b+f^--!] . . . „bj^„Jj+<^]^ . . . qbj\j[j+cr] 



Og' 



(74) 
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The functions F have the same structure as in (|52D , however, they depend hnearly on 
a part of highest-order derivatives. Here the supposition of the rank for the matrix 

- A={Aj,^,a,\,,G') (75) 



is equivalent to the same supposition for the matrix 



(76) 



Qqb^\.J\^+J] 

Let this rank be equal to [Aj/i]. In this case there exists a final division of indices, 

Oil J (ai| J, Oil j+i) with [oii j] = [Aj/i] , 

such that the equations -Fa 7/, = can be solved with respect to the derivatives 
and we obtain, instead of the two first lines of (0), the following expressions: 

qai\.j[i+J] _ i^a.i\j ^_ _ _ qbj\.j[j+J-i] ^ . . . qbj\Ai+'^-'^] ^ . . . qbj\.j+i[j+J]'^ ^ 

^a,|,[i+(T] _ ^a,|„ ^ . . [i+fT_~l]^ _ _ _ ^b^\„^\j+a]^ _ _ _ gf'j|j[i+<^] j ^ 

2, j = 0, 1, I , a = —I, J — 1, — / <(T_<cr, cr+l<cr+<J— 1. 

Now, let us put together the first two entries of cp""^^'' as ■ ■ ■ g''jkb+°'^^]j — / < a < J and 
remind that for a = J no corresponding cr+ occurs. Furthermore, let us replace the last 
entry of v?"'!" as follows: ■ ■ ■ g^^vb+'^l -> ■ ■ ■ q^j\A3+^] . . . qhj\j+i[j+<r] ^ -/ < < J - 1, then 
we get the missing contribution to (T+ for the case under consideration. So, we end up 
exactly with Eq. ( |5^ ) and the Theorem 2 is proved. 

If the rank is less than [^j/i] then the above procedure is applied to the functions 
Faj/-^- Doing that we lower the number of the equations that are not yet reduced to the 
canonical form (the equations of the type -Faj/i = ). Remark that such a diminution 
does not happen at the first stage if = . At a certain stage the procedure does 

not lower the number of the above mentioned equations. This can happen when the rank 
of the matrix of the type (|76| ) is maximal, i.e. is equal to the number of the functions of 
the type Faj/^ ■ In such a case we may reduce them to the canonical form as was said 
above. This can also happen when we do not obtain the functions of the type -Faj/j in the 
reduction procedure. That means that already at the previous step the set is reduced to 
the case Af = K — 1, i.e. the possibility of the reduction to the canonical form is proved. 

Finally we stress that the reduction procedure is formulated for sets of equations of the 
type ( ]5^ ) (the ELE are a particular case of such sets). The procedure holds true under 
certain suppositions of ranks. These suppositions demand various Jacobi matrices of the 
type dFs/dq""^^^ to have constant ranks in the vicinity of the consideration point. Here 

= are equations obtained to a given stage of the procedure and g"''' are highest-order 
derivatives in these equations. It is important to realize that proving the equivalence (|56| ) 
we prove at the same time the locality of the operators U and V from (0). In fact, the 
latter proof is provided by the applicability of the Lemmas from the Appendix. 
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5 Gauge identities and action symmetries 



It was demonstrated above that in the general case of singular ELE the number of the 
equations in the canonical form is less than the number of the equations in the initial set 
of the differential equations. This reduction is related to the fact that in the canonical 
form we retain the independent equations only, whereas the initial equations may be 
dependent. The dependence of the equations in the initial set may be treated as the 
existence of some identities between the initial equations. The identities between the 
ELE imply the existence of gauge transformations of the corresponding action. Below we 
discuss this interrelationship in detail. 

First, we introduce some relevant definitions: The relation of the form 

R''Fa = 0, (77) 

where R"^ are some LO, and Fa (g^'') are some LP, is called the identity between the 
equations Fa (g^'^) = 0. The identity sign = means that the left hand side of (^) is zero 
for any arguments gt'^ . 

Any set R = (^R"^^ of LO that obeys the relation ([77|) is called the generator of an 



identity. Whenever R is a generator than nR with some LO n is a generator as well. Any 
linear combination n*Rj of some generators Rj with operator coefficients n* is a generator. 

A generator R will be called nontrivial if the relation^ nR = O (F) can only be provided 
by a LO n of the form h = O (F) . 

A set of generators Rj will be called independent if the relation ra^Rj = O (F) can only 
be provided by ra* of the form = O [F) . Identities generated by independent generators 
will be called independent. 

Note that for any set of LP Fa , there always exist trivial generators. Namely, the 

generators Rtriv = (^R-triv^ — ^ (F) of the form 

R-^.^ = J2 Ff Im^'^I"'^ , = -u^'^"' , (78) 

k,l 

with arbitrary antisymmetric LP li^'^l'^' obviously lead to the identities (ffTl). These iden- 
tities are not, however, connected to the mutual dependence of the functions Fa . 

An independent set of generators R^ is complete whenever any generator R can be 
represented in the form R = A Rg + Rtriv with some LO A . Any two complete sets of 
independent generators Rg and R^ are related as R^ = U^'lig/ + Rtriv , where U is an 
invert ible LO. 

Supposing now that Fa in Eq. (|7^) are functional derivatives of an action. Fa = SS/Sq"", 
such that Fa = are ELE. Let the functions Fa obey all the necessary suppositions of 
ranks such that ELE can be reduced to the canonical form (|5^). Let us write here this 
canonical form as follows^ 

p = q^[i^] _ ^" (. . . . . . = Q , a = (tt, c/) , (79) 



^^By 6(F) we denote LO of the form (|) with aU the LF u^^ = 0{F) , where 

0(F)|^^o = 0. 

^Here, we do not distinguish possible different proper orders of the coordinates. 
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where are gauge coordinates. Moreover, according to the Theorem 2, there exists the 
equivalence 



= U^F, , = K'F, , UlV,^ = 51 , (80) 



where U and V are LO. Now we may consider the identity (|77D as an equation for finding 
the general form for the generator R . Using (RO) we transform this problem to the one 

" a — 

for finding the operators ^ , 

r^.^O, R'' = l\\ (81) 

Using the explicit form ( pO] ) of the functions Fa , we get ^ = ^ j ; = («; S') ; 

where ^ obey the equation 

fr^O, (82) 



and is a set of arbitrary LO. Since the functions / have the canonical form (pO]), any 
solution of the equation (|8^) is presented by trivial generators of the form 



k,l 



where y}-°''^^°' are arbitrary antisymmetric LF. To demonstrate that we present the gen- 
erators ^ as ^ = J2k=o^°' d'^/df' , where are some LF. Then, in the equation 
(^), we pass from the variables g^H, k,l = 0,1,... to ones g"!'^"], g^H, 

ka = 0, 1, — 1, I = 0,1,... . Such a variable change is not singular. In terms of 
the new variables, the equation (^) reads 



K 



Its general solution is well known 

£-afe \ ^ fV] la'\ka la'\ka ka\la' 

S — / , J a'"- ' " — U 

I 

Now we can write the general solution of the equation ( PT] ) as: 



k,l ^ ^ 



, m'^I'^'^ = (84) 



Let h = {a', g') , a = (a, ^f) in (|^) . Then -u's'l^", li'^'I'^s = —ykaW g^,^^ yig'\kg arbitrary 
LF (e.g., they can be selected to be zero). Indeed, the functions v}^ ^^'^ and v}^ do not 
enter the exjjressions for the generators ^ . Besides, terms with v}°''^^^ affect only the 
generators ^ , which are arbitrary by the construction. Respectively, the general solution 
of the equation ( |77[ ) reads: 

R = f R, + Rtnv , R3 = {Rl = 5lV,- = Vg) , (85) 
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and 



E 

k,l 



U 



lb\kd 



i;") is 



where T^^\^°' = — J'^'^l'^ are some LF. The set of the generators = {Rg 
complete and it is presented by LO. Moreover, these generators are independent. Indeed 
multiplying the equation h^Rg = O (F) from the right by ?7^, we get: h^Sl = 6 [F) 
n9 = 0(F). 

Thus, there exist the following nontrivial identities between the ELE: 

6S_ 



= 0, g 



with generators that are LO. These identities are called the gauge identities. As is 
well known (see for example, 0, §), the existence of the gauge identities (|86D implies the 
existence of infinitesimal gauge transformations of the form 



q" + Sq" ,Sq''={lf] e 



where r parameters = e^(t) are arbitrary functions of time t. Remark that R"^ are LO 
as well. 

Thus, it was demonstrated that for theories that obey appropriate suppositions of 
the ranks there exists a constructive procedure of revealing the gauge generators. For 
such theories all the generators are LO. The number of the independent generators and, 
therefore, the number of the independent gauge transformations is equal to the number 
of the gauge coordinates in the ELE. 

As a simple mechanical example, consider the action of the formQ 



S 



Ldt. 



X' 



The corresponding ELE are: 

Fi = X — y + ax = 0, F2 = x — {I + a)y 
where F2 = is a constraint. The generalized Hessian reads: 



M 



^ = 1 



dydx 



-1 



dxdy 



= -1 

a + 1 



^9) 



(90) 



Let a 7^ , M 7^ 0. In such a nonsingular case the reduction procedure looks as follows: 
By the help of the consistency condition F2 = =^ x = {1 + a)y , we eliminate x from 
the first ELE. Thus, we get an equivalent set, which has the canonical form. 



y 



-X , X 



[l + a)y. (91) 



Another canonical form 



X 



-{l + a)x, y 



1 + ay^x, (92) 

''At a ^ we have a finite-dimensional analog of the Proca action, and at a = we have the analog 
of the Maxwell action. 
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we obtain eliminating y from the equation Fi = by the help of the consistency condition 
= ^ y = x/ (1 + a) . 

Let a = . The case is singular, M = , and rank of the Hessian matrix is equal to 1 . 
One can easily see that the equivalence 

\ _ y( X-y\ y_f d/dt 1 \ ^1 _ / 1 

F2 )~ \ )' \ 1 Q )' ~\l ~d/dt 
holds true. Then the canonical form of the ELE reads x = y and there is a gauge identity 

R^Fa = 0, = 1, = -d/dt . 

The operators transposed to are (^R^^ = ^(-^^) ~ ("^'^) ~ ' "^^^^5 at a = 0, 

the action (|8^) is invariant under the gauge transformation x ^ x + e , y y + e . In the 
case under consideration, the ELE have two canonical forms: x = y and y = x . 



6 Concluding remarks 

We have formulated the reduction procedure which allows one to transform the ELE to 
the canonical form as well as to establish possible gauge identities between the equations. 
The latter part of the procedure can be considered as a constructive way of finding all 
the gauge generators within the Lagrangian formulation. At the same time, it is proven 
that, for local theories, all the gauge generators are local in time operators. The canonical 
form of the ELE reveals their hidden structure, in particular, it presents the spectrum of 
possible initial data, and it allows one to separate coordinates into nongauge and gauge 
ones. One ought also to remark that the reduction procedure can be, in particular, treated 
as a procedure of finding constraints in the Lagrangian formulation. 

In that respect one can compare the reduction procedure with the well-known Dirac 
procedure in the Hamiltonian formulation of constrained systems 0, ^. Recall that 
the Dirac procedure is applicable to the Hamilton equations with primary constraints, 
namely to equations of the form 

F{j^^r]) =7]- [t], H^^^ } = , <l>(^) {v) = 0, H^^^ = H [r]) + A$(^) (r/) . (93) 

Here r] = {q"',Pa) are phase-space variables; ^^^^ {r]) = are primary constraints, A's are 
Lagrange multipliers to the primary constraints, and H^^^ is the total Hamiltonian . Via 
{■, ■} the Poisson bracket is denoted. The aim of the procedure is to eliminate as many as 
possible A's from the equations, to find all the constraints in the theory. The procedure is 
based on the consistency conditions <i>^^^ = 0. Using the equations F (r/, r/) = 0, we may 
transform any consistency condition to the following form: 

<i.« = {$«,i7(i)} = 0. 

From these equations one can define some A's as functions of 1] and reveal some new 
constraints. Then the procedure has to be applied to the latter constraints and so on. 

The equations (0) present a particular case of differential equations considered in 
the present article (indeed, these equations are ELE for a Hamiltonian action). Thus, 
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our reduction procedure may be applied to these equations. Namely, first one has to 
consider the equations Fa = 0, = and select independent w.r.t. f] equations. Since 
equations of the primary constraints are independent by construction, we pass to the next 
step and solve the constraint equations ^^^^ = with respect to a part of the variables r], 
as = 

Vi- {1I2) = 0- Then we exclude 77^ and r/^ from the equations F = 0. 
Thus, we get F = — > FA{ri2,ri2) = 0. Then one has to select independent w.r.t. r/2 
functions Fa^^ ■ At the same time one finds new constraints Fa^^ (^2) = so on (see 
the Subsect. "First order equations"). 

We see that the Dirac procedure differs from our reduction procedure. Indeed, as was 
mentioned above, in the Dirac procedure one excludes all the derivatives r/ by the help of 
the equations F = from the consistency conditions = 0. Thus, one gets equations 
for the Lagrange multipliers and new constraints. Besides, one of the aim of the Dirac 
procedure is to maintain the canonical Hamiltonian structure of the equations F = 0. The 
possibility of the Dirac reduction is due to the specific structure of the equations (|93D. 
Namely, here the consistency conditions never involve A and rajik dF/di] = [F] = [77]. 



Besides, one ought to mention the work where it was proposed an alternative (to 



the Dirac procedure) way of reducing the equations of motion for theories with actions 
of the form S = J [ipA (v) V"^ — V (r^)] dt . One can verify that, in fact, the procedure of 
that work, in a part (the procedure does not reveal the gauge identities), is similar to our 
reduction procedure in the case of the first order equations (see Sect. IV). 

However, the reduction procedure proposed in the present article is formulated for a 
wider class of Lagrangian systems (differential equations). It does not need the intro- 
duction of new variables such as momenta and Lagrange multipliers, and it is defined 
in the framework of the initial Lagrangian formulation. Moreover, its aim is twofold: to 
reduce ELE to their canonical form and to reveal the gauge identities between the ELE 
equations. 

The consideration in the present article is restricted by finite-dimensional systems. Its 
application to field theories (theories with infinite number degrees of freedom) demands 
additional study. We hope to present the corresponding formulation in futures publi- 
cations. However, in simple cases, one can apply the present reduction procedure with 
some natural modifications in the infinite-dimensional case. Consider the Maxwell action 
S = — (1/4) f J^^^J^^^^dx , J-'^y = d^Ay — d^A^ as a common example of a singular field 
theory. The ELE read: 

F' = —- = d.r" = A' + diA^ - AA' + d,^ = Q , (94) 
oAi 

F° = -^ = a.^-o = (^ + AAO =0, ^ = afcA^ (95) 

The equation Fq = is a constraint. Following the reduction procedure, we have to 
consider the set Fj = , Fq = . The Jacobi matrix dF^/dA^, has the constant rank 
3. We can, for example, select the equations (0) as independent with respect to the 
derivatives A*. The equation Fq = Ois their consequence. No more constraints appear. 
Now we exclude A^ and A^ from (p^ by the help of (P^). That creates the equivalence 
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Now we discover that the functions are dependent, d^F^ = 0. In our terms that reads, 
for example, as the following equivalence 



( F^\ I 1 \ / F 

FM = 1 

\F^ I \ -d^'d, -d^'d2 1 / V 

The equations F^ = 0, = 0, F° = present one of the canonical forms of the Maxwell 
equations. The identity that follows from the presence of the zero in the right column of 
(p7|) reads as S^-F'^ = in terms of the initial functions F'^ and implies the invariance of 
the Maxwell action under gradient gauge transformations. 
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7 Appendix 

Here we present three Lemmas which are used in the reduction procedure to justify 
equivalence of equations and LF. In this respect it is useful to recall here the relevant 
definitions from the Sect. II . Namely: 

Two sets of equations. Fa (g''^) = and fa {q^''^) = are equivalent F = <^=^ / = 
whenever they have the same set of solutions. If two sets of LF Fa (q'''^) and Xa ' 
[F] = [x] , are related by some LO U and as F = Ux , X = ? UV = 1 , then we 
call such LF equivalent and denote this fact as: F ~ x . In this case the corresponding 
equations are strongly equivalent. 



Lemma 1 

Let a set of equations 

(x, y =0, Fa{x,y)=0, x = (x^) , y = (y") , (98) 

be given, where $ are some LF. And let det dFa/dy''\^^y^ ^ 0, where the consideration 
point {xo,yo) is on shell. Then: 

a) The equations Fa {x,y) = can be solved w.r.t. y as: y"" = (p"'{x) , where ip°'{x) are 
some single-valued functions of x in the vicinity of the point xq . In other words, there is 
the equivalence 

Faix,y)r~.y'^-cp''{x), (99) 

which implies the strong equivalence between the equations F^ (x, y) = and y"" = (p°'{x) . 

b) The following equivalence between the LF holds true: 



Fa{x,y) I '\y--^''{x) I' - "^Ml,w=^m • (100) 
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The first statement is, in fact, tlie well-known implicit function theorem |13|. Taking 
into account (pOf), we have: Fa {x,y) = Uab iy"' — '{'"'{x)) , det ^^Lgyg 7^ 0. On the other 
side one can write $^ = $^ + V^a [y"" — '^"'{x)] , where V^a is a LO. Thus, 



F 



U 



, u 



1 V 

M 



1 -1/ 

' ' 



(101) 



and the equivalence (|10CI|) is justified. 

Lemma 2 

Let a set of equations 

F4 (g, = , g = {q") , z = (2;') , A = 1, m, a = 1, n , i = 1, / , 

be given. And let the Jacobi matrix dFA/dq"' have a constant rank in a vicinity Do of the 
consideration point (go , -Zo), which is on shell {Fa (go , -^o) = 0), 



rank 



BFa 



r . 



(102) 



Then there exists an equivalence 

yf^ — ift^^X, Z) 



, q- = {x\yn, A=(/i,G), [/i]=r. 



(103) 



We begin the proof with the remark that due to ( |102| ), there exists a division of the 
indices A = (/x, G), a = (/x , (7) , [/x] = r , g" = (x^, j/'^) , such that 



det ^ 



7^0 



90,^0 



Then by virtue of the Lemma 1 we can write 



u^,u^, r = y'' - V^ix, z), det wj^^^^^ ^ . 



(104) 



(105) 



Let us present the functions Fq in the form Fq {x, y, z) = Qq {x, z) + 11^/^/'^ {x, y, z) , 



where VLq (x, z) = Fg\ 



y=ip{z,x) ' 



such that ^Ig {xq, zq) = . Then 



UabXb 5 Xb 



no 



u 



u 

n 1 



detf/|,„,„7^0. (106) 



In virtue of (|102|) and ([IO6D 



rank 



9g" 



(107) 



q,z&Do 



The Jacobi matrix 9xy^/9g'*has the following structure: 



Sg'^ 9(?/'',a;9) 



-d^^jdx^ 
dVtc/dxa 
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Therefore, 



rank 



one 

dx3 



one 

dx3 



x,z^Do 



Eq. (|108| ), together with the relation {^o, ^o) = 0, imphes 



Finally, we may write 

Fa = UabXb 5 Xb 



Thus, the equivalence (p.03|) is justified 



f {x, y, z) 
^g{z) 



Lemma 3 



det f/l 

1 90,^0 ' 



(108) 



(109) 



As a consequence of the Lemma 2 the following Lemma holds true: 
Let a set of equations 

Fa {q") = , A = 1, m, a = 1, n , 

be given. And let the Jacobi matrix dFA/dq"" have a constant rank in a vicinity Dq of the 
consideration point go which is on shell (F (go) = 0), 



rank 



r . 



qeDo 



Then there exists an equivalence 
Fa ~ Fa = 



Og 



, A={fi,G),OG = OWG, [fi]=r 



(110) 



The proof of this Lemma follows the one of the Lemma 2 if one selects there z = zq. 
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